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Abstract 

We compare first order (refractive) ionospheric effects seen by the Murchison Widefield Array (MWA) with 
the ionosphere as inferred from Global Positioning System (GPS) data. The first order ionosphere manifests 
itself as a bulk position shift of the observed sources across an MWA field of view. These effects can be 
computed from global ionosphere maps provided by GPS analysis centres, namely the Center for Orbit 
Determination in Europe (CODE), using data from globally distributed GPS receivers. However, for the 
more accurate local ionosphere estimates required for precision radio astronomy applications, data from local 
GPS networks needs to be incorporated into ionospheric modelling. For GPS observations, the ionospheric 
parameters are biased by GPS receiver instrument delays, among other effects, also known as receiver 
Differential Code Biases (DCBs). The receiver DCBs need to be estimated for any non-CODE GPS station 
used for ionosphere modelling, a requirement for establishing dense GPS networks in arbitrary locations in 
the vicinity of the MWA. In this work, single GPS station-based ionospheric modelling is performed at a time 
resolution of 10 minutes. Also the receiver DCBs are estimated for selected Geoscience Australia (GA) GPS 
receivers, located at Murchison Radio Observatory (MROl), Yarragadee (YAR3), Mount Magnet (MTMA) 
and Wiluna (WILU). The ionospheric gradients estimated from GPS are compared with the ionospheric 
gradients inferred from radio source position shifts observed with the MWA. The ionospheric gradients at all 
the GPS stations show a correlation with the gradients observed with the MWA. The ionosphere estimates 
obtained using GPS measurements show promise in terms of providing calibration information for the MWA. 

Keywords: atmospheric effects - techniques: interferometrici 
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1 INTRODUCTION 

The Earth’s ionosphere, being a dispersive medium at 
radio wavelengths, causes a change in the propagation 
velocity of radio waves, among other effects. The 
ionosphere is directly dependent on the solar activity 
through the high energy far ultraviolet and X-rays. 
The ionosphere further varies depending on various 
transportation and depletion processes, namely, due to 
the influence of tides and atmospheric (gravity) waves, 
solar winds, and vertical transport through eddy’s 
diffusion and the geomagnetic disturbances (Zolesi 
& Cander, 2014). Considering the ionosphere to be 
highly variable and uncertain, it becomes important to 
monitor the ionosphere on a regular basis for a number 
of applications. 

Early ground-based ionospheric sensors like ionoson- 
des were effective for understanding the bottomside 
ionosphere. Ground based Coherent and Incoherent 
radars operating at High Frequency (HF, 3-30 MHz), 
Very High Frequency (VHF, 30-300 MHz) and Ultra 
High Frequency (UHF, 300-3000 MHz), are able to 
probe the middle and upper ionosphere. Example of 
coherent radars include the Super Dual Auroral Radar 
Network (SuperDARN) consisting of around 34 HF 
radars located at mid-latitudes and extending into 
the polar regions. The radars are pointed towards 
the North and the South poles to study ionospheric 
convection (Zolesi & Cander, 2014). Incoherent radars 
include the Jicamarca Radio Observatory located 
along the geomagnetic equator in Lima, Peru, (49.92 
MHz), the Arecibo dish in Puerto Rico (430 MHz), and 
EISCAT (European Incoherent SCATter), northern 
Scandinavia, operating at UHF (931 and 500 MHz) and 
VHF (224 MHz) (Zolesi & Cander, 2014). To access 
information about the topside ionosphere, satellite 
based topside sounders have also been used (Hun- 
sucker, 2013). With the advent of radio communication 
satellites, the signals from such satellite systems were 
used to obtain spatial and temporal information on the 
ionosphere (Leitinger et al., 1984). 

The Global Positioning System (GPS) was designed 
by the Department of Defence, U.S., in the early 1970s 
to fulfil U.S. military requirements (FRP, 2001). It has 
since been used in various civilian applications, from 
navigation to precise geodetic positioning. The promise 
of GPS to operate in all weather conditions, 24x7, in 
addition to its multi-frequency transmission, has made 
it a useful tool to monitor ionospheric parameters 
(El-Rabbany, 2006). In comparison to ground-based 
ionospheric sensors like radars and ionosondes, the 
satellite based system like GPS can provide continuous 
near real-time global coverage of the ionosphere. 


For precise positioning applications using GPS, 
new control points can be established by applying 
constraints using positions from the established control 
points. Global GPS data processing centres like the 
International GNSS (Global Navigation Satellite Sys¬ 
tem) Service (IGS), with international multi-agency 
members, provide support for such global geodetic 
activities. In addition to the positions of the global 
network of GPS/GNSS stations, a number of products 
such as precise satellite ephemerides, satellite clock pa¬ 
rameters, Earth rotation parameters, global ionosphere 
maps, and zenith tropospheric path delays are routinely 
generated (Beutler et al., 1999). Global ionospheric 
maps are also generated by various IGS analysis 
centres, namely, the Center for Orbit Determination 
in Europe (CODE), the Astronomical Institute at the 
University of Bern (AIUB), Switzerland (Rothacher 
et al., 1997), and the Jet Propulsion Laboratory (JPL), 
California Institute of Technology, U.S (Mannucci 
et al., 1998; Komjathy et al., 2005), among others. 
The temporal and spatial resolution of global maps is 
generally of 2 hours and 5°/2.5° in longitude/latitude, 
respectively. 

Real-time GPS positioning accuracy can be im¬ 
proved by providing ionospheric and other corrections 
directly to the user. This is the role of Space Based 
Augmentation Systems (SBAS) such as the Wide Area 
Augmentation System (WAAS) in the US (Parkinson 
& Spilker, 1996). Vertical ionosphere corrections are 
generated at temporal and spatial resolution of three 
minutes and 5°/5° in latitude/longitude, respectively 1 . 

The Global GPS data are also processed at the 
MIT Haystack Observatory by the MAPGPS software 
package in order to generate global TEC maps (see 
Rideout & Coster, 2006). The TEC maps are generated 
with a greater temporal and spatial resolution of 10 
minutes and 1°/1° in longitude/latitude, respectively 
and are distributed through an open source, web 
based system 2 . To achieve even finer spatial resolution, 
regional ionosphere modelling must be performed. 
Regional ionosphere data centres like the Royal 
Observatory of Belgium (ROB) are able to model 
the ionosphere over small spatial areas with higher 
sampling of ground GPS stations with a temporal 
and spatial resolution of 15 minutes and 0.5°/0.5° in 
longitude/latitude, respectively (Chevalier et al., 2013). 

Since the early days of radio astronomy, the iono¬ 
sphere has been found to have a profound effect on 
astrometric observations. In a study of discrete radio 
sources, fluctuations in the signal from Cygnus A 
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at observing frequency of 68 MHz were reported by 
Hey et al. (1946). Further investigation confirmed 
this to be a result of ionospheric structures. Hewish 
(1951) and Booker (1958), among others, used this 
information to quantify ionospheric fluctuations and 
present some insights into the nature and behaviour of 
the ionosphere. Radio astronomy can be used to extract 
information on the ionosphere, however an external 
source for deriving information on the ionosphere is 
needed to calibrate astrometric observations. 

Ros et al. (2000) evaluated the quality of GPS-based 
ionosphere corrections for Very Long Baseline Interfer¬ 
ometry (VLBI) at 8.4 and 2.3 GHz. In his study, Ros 
et al. corrected the ionosphere for continental (~200 to 
~700 km) and inter-continental (~7800 to ~8300 km) 
baseline lengths and concluded that GPS maps of Total 
Electron Content (TEC) can usefully contribute to 
VLBI astrometric analysis. Erickson et al. (2001) made 
use of the ionospheric corrections generated from an 
experimental set-up of four GPS receivers at the Very 
Large Array (VLA) site to correct for the ionospheric 
effect. Erickson et al. found that large scale structures 
(>1000 km) could be resolved observing at frequencies 
of 322 and 333 MHz, whereas small scale fluctuations 
(<100 km) could not be seen using a global model; 
it was noted that global ionosphere models perform 
averaging over the ionosphere and hence lose their 
capacity to monitor small scale ionospheric changes. 
Erickson et al. argues that a dense GPS network is 
required to correct for small scale fluctuations in the 
ionosphere. 

In work by Sotomayor-Beltran et al. (2013), CODE 
and ROB maps were used to compute the Rotation 
Measure due to the ionosphere for polarised sources 
that were observed by LOFAR (the Low-Frequency 
Array for radio astronomy). However, the CODE maps 
presented in Sotomayor-Beltran et al., (Figure 2) show 
an error in the location of the equatorial anomaly. This 
error is further investigated and discussed in Appendix 
A of this paper. This error is also discussed by Herne 
et al. (submitted, A & A). 

GPS-based estimation of the ionosphere has also 
been carried out by Herne et al. (2013). High-fidelity 
GPS systems were specially deployed by the US Air 
Force for measurement of ionospheric TEC and scintil¬ 
lation indices at the location of the Murchison Radio 
Observatory (MRO) and Australian Space Academy 
campus (Meckering, Western Australia). Two periods 
of time corresponding to low (2008-2009) and high 
(2012-2013) ionospheric activity were studied. During 
2008-2009, the F10.7 Index varied between 65 to 82 
(10 ^ 22 Wm~ 2 Hz^ 1 ), the monthly variation was found 
to be insignificant for most of the period. The Ap index 


reached a maximum of 35 during this period. However, 
during high solar activity (2012-2013), F10.7 led be¬ 
tween 85 to 190, the monthly variation was significant 
throughout the period. Ap index reached a maximum 
of 90. The ionosphere at the MRO was expected to 
exhibit very low levels of scintillation (Kennewell et al., 
2005). For Kp=9, scintillation at MRO reached about 
0.2. In Herne et al., during 2012-2013, scintillation was 
found to reached a maximum of about 0.3. 

The Murchison Widefield Array (MWA) is a recently 
operational low-frequency radio telescope located 
within the MRO in Western Australia. The instrument 
exemplifies exceptional wide-held imaging capability 
at low frequencies (see Lonsdale et al. (2009) and 
Tingay et al. (2013) for a technical overview of the 
instrument). These capabilities make it ideal for 
a range of science investigations with an emphasis 
on: detection of the 21cm line from the Epoch of 
Reionoisation; large-scale surveys; searches for radio 
transients; continuum surveys; and Solar, Heliosphere 
and Ionospheric studies (Bowman et al., 2013). In the 
recent survey by Hurley-Walker et al. (2014), over 
14,000 radio sources were detected in just four nights 
of observing. 

Since a fundamental observable of an interferom¬ 
eter is the phase difference between elements, such 
instruments are unable to measure the Total Electron 
Content (TEC) towards a particular source directly 
(though it may be inferred indirectly from polarisation 
measurements). However, for a single source located 
at infinity, the radiation reaching each element of an 
interferometer will have passed through a different 
part of the ionosphere providing information on the 
differential ionosphere or ionospheric gradients. 

Excellent summaries of how the ionosphere affects 
radio interferometers are given in Lonsdale (2004) and 
Wijnholds et al. (2010). In the specific case of the 
MWA, the instrument has a very wide field of view, but 
relatively short baselines. Thus, while the effects of the 
ionosphere may change across the instruments’ field 
of view, the complex structure of the ionosphere is on 
larger spatial scales than the dimensions of the array, 
and the wavefront arriving at the array from a source 
will be largely coherent. However more complicated 
effects may be seen. Unless explicitly corrected in 
the calibration process, a TEC gradient across the 
array will manifest itself as a shift in source position; 
the incoming radiation gets refracted by the ionosphere. 

The MWA has recently commenced operations and 
ionospheric effects are routinely detected. For example, 
observations have revealed variations in the Rotation 
Measure of polarised point sources and the diffuse 


PASA (2015) 

doi: 10.1017/pas. 2015.xxx 


4 


Arora et al. 


galactic background (Lenc, private communication). 
Loi et al. (2015a) demonstrate the utility of the MWA 
as a powerful imager for studying high-altitude irregu¬ 
larities, revealing a population of field-aligned density 
ducts that appear regularly over the observatory. The 
spatial distributions of celestial source refractive offsets 
over the field of view, combined with a novel parallax 
technique for altitude measurement, enable the 3D 
characterisation of ionospheric structures at high 
temporal cadence. Travelling ionospheric disturbances 
(TIDs), sinusoidal perturbations with wavelengths of 
100-1000 km and periods of several tens of minutes 
to an hour, are also often observed in MWA data. 
A technique for the spatio-temporal power spectrum 
analysis of ionospheric gradients, presented by Loi 
et al. (2015b), enables characteristic wavelengths and 
periods of fluctuation to be measured. A quantitative 
study of the statistical properties of ionosphere-induced 
position and amplitude variations of sources in MWA 
images (Loi et al., submitted. 2015) has yielded infor¬ 
mation about characteristic density gradients and the 
diffractive scale in the ionosphere. This illustrates the 
capability of the MWA to probe ionospheric structure 
in great detail. 

All previous attempts at GPS-based ionosphere 
calibration for radio astronomy made use of ionosphere 
information from TEC maps which have a temporal 
resolution of 2 hours. Ionospheric structures such as 
TIDs have a period of several tens of minutes, hence 
it is important to generate ionospheric information 
with higher resolution. We present here a method 
where publicly available data from GPS stations are 
used to generate location-specific ionosphere models at 
10 minutes intervals. This requires the calibration of 
GPS receiver/transmitter instrumental biases in order 
to accurately model the ionosphere, also known as 
Differential Code Bias (DCB). 

We perform a simultaneous fit for both GPS DCBs 
and ionospheric parameters. While the Bernese 
software provides the GPS DCBs and ionospheric pa¬ 
rameters (Beutler et al., 2007), our work allows much 
higher time resolution of the ionospheric parameters 
than is commonly obtained, while simultaneously 
solving for the DCBs. 

This paper is organised as follows. In section 2, 
the methodology concerning the GPS-based ionospheric 
modelling is fully described. In this section, the basic 
GPS observation equations, estimating the Vertical To¬ 
tal Electron Content ( VTEC ) are formulated. In par¬ 
ticular, we show how the VTEC and DCBs are deter¬ 
mined simultaneously through Kalman filtering. Section 
3 briefly describes the methodology adopted by CODE 


to generate ionosphere maps. Section 4 presents the pro¬ 
cedure to model the ionosphere seen by the MWA in the 
form of position shifts. The results are presented and 
discussed in section 5. Results from the IGS analysis 
centre, CODE, serve as reference to validate our VTEC 
results. In particular, we compare our estimated receiver 
Differential Code Biases (DCBs), the inter-frequency 
biases on code observables, with those determined by 
the Bernese processing software employed by CODE. 
The ionosphere observed by the MWA and the GPS are 
analysed and the agreement discussed. Finally, we make 
some concluding remarks in section 6. 

2 ESTIMATION OF THE IONOSPHERE 

USING GPS OBSERVATIONS 

A GPS constellation consists of up to 32 operational 
satellites placed in 6 orbital planes at an altitude 
of 20,200 km above the Earth’s surface. With this 
constellation geometry and an orbital period of about 
12 hours, 4 to 10 GPS satellites are visible anywhere 
in the world at any given time. The earlier GPS 
satellites (Block IIA (2nd generation, “Advanced”) 
and Block HR (“Replenishment”) transmitted on two 
carrier frequencies (LI and L2), each encoded with 
one or two digital codes, and navigation messages (El- 
Rabbany, 2006). The BLOCK IIR(M) (“Modernized”) 
satellites have an additional civil signal on L2 (L2C) 
(Hofmann-Wellenhof et al., 1993). The BLOCK IIF 
(“Follow-on”) satellites transmit an additional third 
frequency, L5, which has higher transmitted power 
and greater bandwidth, to support high-performance 
applications (El-Rabbany, 2006). 

The signals transmitted by the GPS satellites form 
the observables, namely the phase (of the carrier 
frequency) and code (digital code) measurements. A 
phase measurement is the number of cycles at the cor¬ 
responding carrier frequency between the satellite and 
the receiver. The phase delay between the receiver and 
the satellite is obtained by multiplying the number of 
phase cycles with the wavelength of the corresponding 
carrier. The phase measurements are biased by an 
unknown number of phase cycles, in addition to other 
errors. When a GPS receiver is switched on or tracks 
a newly risen satellite, it cannot determine the total 
number of complete cycles between the receiver and the 
satellite. Hence the initial number of complete cycles 
remains ambiguous, known as phase ambiguity bias. 

The GPS code measurements have two types of code 
observables, namely the C/A-code (Coarse/Acquisition 
code, modulated only on the LI carrier, denoted as Cl) 
and P-code (Precise code, modulated on both LI and 
L2 carriers, denoted as PI and P2, respectively). The 
code modulation is different for each GPS satellite, with 
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code signals sometimes also referred to as PRN (Pseudo 
Random Noise). The C/A-code measurement is less pre¬ 
cise than the P-code, since the bit rate of C/A-code is 
10 times lower than P-code (Langley, 1993). The GPS 
receiver generates replicas of the transmitted code sig¬ 
nals. By comparing the code signal to its replica, the sig¬ 
nal travel time is obtained (Hofmann-Wellenhof et al., 
1993). For a more detailed description, one can refer to 
Hofmann-Wellenhof et al. (1993), Langley (1993) and 
El-Rabbany (2006). This work uses two GPS carrier 
frequencies, namely /i = 1575.42 MHz (for carrier LI) 
and f 2 = 1227.60 MHz (for L2). 


2.1 The GPS observation equation 

The terms contributing to the GPS phase and code ob¬ 
servables can be given as follows (Teunissen & Kleus- 
berg, 1998): 


C* = [c • (S r , 2 i - J‘ 21 ) + (Ai M‘ tl - A 2 M“ 2 )] . (6) 

The slant ionospheric delay, •, is related to the 
Slant Total Electron Content ( STEC ) as follows 
(Hofmann-Wellenhof et ah, 1993): 


, 40.3 STEC 

r,j f 2 

J 3 


( 7 ) 


As shown in Figure 1, the STEC can be mapped 
to the Vertical Total Electron Content, VTEC , using 
the obliquity factor F s , also known as the ionosphere 
mapping function. F s is a function of zenith angle at 
the Ionosphere Pierce Point (IPP), z', given as follows: 


STEC = VTEC ■ F s 


E(K,i) =Pr~ + c • (S r j - S'j) + A jM’j, (1) 

E{P^)=-pl + ^ j +c-{d r , j -d° j ). (2) 

The parameters used in sections 2.1 and 2.2 are listed 
in Table 1. Here subscripts s,r, and j indicate satellite, 
receiver and GPS frequency number, respectively. 

Af/ ■ are the non-integer ambiguities on the phase ob¬ 
servables which contain the unknown integer ambigui¬ 
ties, V/ j , and the non-integer initial phase offsets for 
the receiver and satellite he., 

M r,j =A r r,i + ^ry(t0)-^(t 0 ). (3) 

The phase ambiguities remain constant for any given 
receiver, frequency, and continuous satellite arc unless 
there is a loss of signal lock. 

Since the effect of the ionospheric delay is a func¬ 
tion of frequency, a frequency difference of phase and 
code observables can be formed which retains the iono¬ 
spheric delay while the geometry-related parameters are 
eliminated. Along with the ionospheric delay (t*), the 
phase and code instrumental delays (S r j,Sj, d r ^,d s ^) 
and phase ambiguities remain. The frequency- 

difference phase and code observation equations for dual 
frequency GPS observables, formed using equations (1) 
and (2), are given as follows: 


F s 


1 


sm z 


11 ' \J 1 — sin 2 z' 

Re 


Re + R-ion 


■ sin(z) 


( 8 ) 


where z' is the zenith angle at the IPP, R e is the mean 
radius of the Earth, considered to be 6371 km and as¬ 
suming a spherical Earth, Hi on is the height at the sub- 
ionospheric point, assumed to be 450 km, and z is the 
zenith angle of the satellite as seen by the receiver. The 
geometry of the model is illustrated in Figure 1. This 
study aims to compare and analyse ionosphere gradi¬ 
ents, the height, Hi on , of 450 km was chosen in order to 
compare the VTEC with CODE analysis center pub¬ 
lished values. 

We can now map the GPS observables to the VTEC 
as follows 


/ $;,2i \ 

V40-3^2i/ 

f-SsU 

V40.3/T21 J 


=-F s VTEC 

= F s VTEC + 


Ct 


40.3 p,2 1 


C- (dr,21 - d* 21 ) 
40.3 /i 2 1 


>( 9 ) 


Here c • (d r , 2 i) and c • (d s 2 i) constitute terms known as 
Differential Code Bias (DCB) for the receiver and the 
satellite, respectively, where 


^r,2l) = ^,l-^,2 = -<21+C? 


( 4 ) 


E(P ^ 2i ) = P/ ;1 - P/, 2 = t S r>21 + C ■ (dr,21 - df 21 ). (5) 

Here, d s 21 and d r , 2 i are the inter-frequency code de¬ 
lays for the satellite and the receiver and the constant 
phase term C* is given as follows: 
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Ml - 72 P2 ~ -72 
J 1 J 2 


(10) 















6 


Arora et al. 


Table 1: GPS model parameters and definitions 


Parameter 

Definition 

Units 

E(-) 

Expectation operator 


cb s . 
r,j 

Phase observables 

(metres) 

P r,J 

Code observables 

(metres) 

Pr 

‘Geometry’ parameters, p* + r* + c • dt* 

(metres) 

Pr 

Receiver-Satellite range 

(metres) 

T S 

1 r 

Tropospheric path length 

(metres) 

c ■ dt s r 

Receiver and satellite clock errors 

(metres) 

b r ,3 

Receiver-Satellite ionospheric delay 

(metres) 

Pi = 1 /fj 

Inverse of frequency square 


fi 

GPS frequencies, j = 1,2 


M s • 

r,j 

Non-integer phase ambiguity 

(cycles) 

N s • 

r,j 

Integer phase ambiguity 

(cycles) 

*Pr,j ) 

Initial receiver phase offset 

(cycles) 


Initial satellite phase offset 

(cycles) 

A, 

wavelength at frequency j 

(metres) 

c-dfi 

Code satellite instrument delay 

(metres) 

c ■ d r j 

Code receiver instrument delay 

(metres) 

c ■ s j 

Phase satellite instrument delay 

(metres) 

c * 

Phase receiver instrument delay 

(metres) 

cb s 
^r, 21 

Frequency-difference of phase observables 

(metres) 

P s 

r r, 21 

Frequency-difference of code observables 

(metres) 

t r,21 

Frequency-difference of ionospheric delay 

(metres) 

c° 

Constant term over each satellite arc 

(metres) 

C ■ d r ,21 

Receiver DCB 

(metres) 

c ■ d s 2 1 

Satellite DCB 

(metres) 

STEC 

Slant Total Electron Content 

(TECU) 

z 

Zenith angle of satellite 

(radians) 

z' 

Zenith angle at IPP 

(radians) 

R e 

Radius of Earth 

(metres) 

Hion 

Height of the ionospheric layer 

(metres) 

F s 

Mapping function 


P21 = Pi ~ P2 

Frequency-difference of 


VTEC 

Vertical Total Electron Content 

(TECU) 

\Pm 

Geomagnetic latitude at IPP 

(radians) 

S 

Sun-fixed longitude at IPP 

(radians) 


Geomagnetic longitude at IPP 

(radians) 

VTECq 

VTEC at receiver location 

(TECU) 

Pm 0 

Latitude at receiver location 

(radians) 

so 

Sun-fixed longitude at receiver location 

(radians) 

S' 

First order derivative 


f" 

Second order derivative 


DCB S 

Receiver DCBs lumped with satellite DCBs(= c • (d r> 21 — 

^?2l)) 

(metres) 


2.2 GPS ionospheric modelling 

The VTEC can be modelled by assuming that the iono¬ 
sphere is concentrated at a single layer at height Hi on as 
illustrated in Figure 1 (Schaer, 1999; Wang et al., 2014). 
The intersection of the GPS receiver-satellite line-of- 
sight with the ionospheric layer is called as IPP. The 


slant TEC is mapped to the vertical TEC by an obliq¬ 
uity factor (equation (8)), which is a function of the 
zenith angle at the IPP, z'. VTEC is modelled as a 
function of geomagnetic latitude, ip m , and Sun fixed 
longitude, s, using the following polynomial function 
(Wang et al., 2014) 
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VTEC(<pm, s ) = VTEC 0 + (ip m - vO/V + 0 - «o)/'s + 

+ (^m - Vm 0 ) 2 f''‘PmVm + (s - S 0 ) 2 /"sS 

+ (^m-¥’m 0 )(s-So)/'W (H) 

The Sun fixed longitude, s, is related to the lo¬ 
cal solar time (LT) as s = A m + LT — 7 r, where A m 
is the geomagnetic longitude at IPP, LT is in radi¬ 
ans, VTECq is the VTEC at the receiver location and 
fs, fipm, f"ss, f"ipm<Pm, f"<Pm.s are the first and 
second order derivatives of VTEC along the Sun fixed 
longitude and latitude, respectively. 

The single layer ionospheric model is a computation 
efficient model to estimate local ionospheric gradients 
at the zenith. However, the ionospheric features, for ex¬ 
ample those along the magnetic field lines, cannot be 
resolved using this model. 



Figure 2.: Flow chart describing the implementation of 
a Kalman Filter. 


m satellites seen by receiver r, a weighted least square 
model is formed using equations (5) and (11), given as 
follows 

E(yi) — Ai Xi 

D(yi) — Qyi (12) 


2.3 VTEC determination through Kalman 
filtering 

The GPS model to estimate unknowns can be formed 
from equation (5). In our method, the DCBs for re¬ 
ceiver and satellite are estimated as a single parameter, 
DCB S = c • (d rt 21 — d s 21 ), since the model presented in 
equation (9) is rank deficient. The parameters, namely 
the receiver and the satellite DCBs, cannot be sepa¬ 
rated from each other, hence cannot be independently 
estimated. A minimum set of parameters, known as the 
5-basis, are chosen which can be lumped with the re¬ 
maining parameters in order to overcome the rank de¬ 
ficiency in the underlying model (Teunissen, 1984). For 


Xi — DCB S VTECq fifim f's f"ss /'Vms] 

where E(-) and D(-) denote the expectation and dis¬ 
persion operators, yi denotes a vector of observables of 
size [ 2 m x 1 ] at time stamp i, Ai is the design matrix 
of size [ 2 m x ( 2 m + 6 )], the unknowns given by Xi are 
of size [( 2 m + 6 ) x 1 ], and Q y i is the stochastic model 
for observables in yi given as 


Qyi — 


40.3/x 2 i / sin 2 (eZ)’ 


(13) 


el is the elevation angle of the satellite, and <To is the 
measurement noise of the observables. 
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To solve the above model, a cut-off for the satellite 
zenith angle of 70° was chosen. The model given in equa¬ 
tion (12) cannot be solved in a single epoch. An iono¬ 
sphere refreshing interval of 20 epochs (1 epoch = 30 
seconds) is chosen over which the ionosphere is assumed 
to remain constant. The phase bias term CJ! constitutes 
an integer phase ambiguity and non-integer receiver and 
satellite phase bias terms. The integer phase ambigui¬ 
ties remain constant for a single continuous satellite arc 
unless there is a loss of signal lock, hence C* is assumed 
to remain constant for a continuous satellite arc. The 
GPS observables were free of multipath and cycle slips. 
The code bias terms, DCB S , are assumed to remain 
constant over a period of 24 hours for any satellite. 

With the above considerations, the unknowns given 
in equation (12) are estimated using the Kalman Fil¬ 
ter approach (Kalman, 1960; Kailath, 1981; Grewal & 
Andrews, 2011). In this work, Kalman filter is used in 
data assimilation mode. A Kalman filter can be de¬ 
scribed as a three step procedure, with initialisation, 
prediction, and measurement update executed at dif¬ 
ferent time steps t, t = 1, •• • ,t max . The implementa¬ 
tion of the Kalman filter is shown in Figure 2, where, 
Xo|o,A q\ q ,Qv, and j/ 0 |o indicate the unknowns, the de¬ 
sign matrix, the Variance Covariance (VC) matrix of 
the measurements, and the measurements, respectively. 
*I> t | t _ 1 is the transition matrix which relates the un¬ 
knowns at the current x t \t-i and previous Xt- i|t_i time 
step with a system noise given by the VC matrix. St, Vt 
are the predicted residuals and Q Vt its VC matrix. K t 
is the Kalman gain which is used to compute the mea¬ 
surement update given by x t \ t and its VC matrix P Xt|t . 
The parameters used in Kalman filter are presented in 
Table 2. 

2.4 GPS DATA PREPARATION 

The data from the three GPS/GNSS stations near¬ 
est to the Murchison Radioastronomy Observatory 
(MRO) were used for this analysis, namely the Geo¬ 
science Australia (GA) stations MROl (Murchison), 
MTMA (Mount Magnet), YAR3 (Yarragadee), and 
WILU (Wiluna) were chosen (Figure 3). A description 
of the selected GPS/GNSS stations is given in Table 3. 
The data for the selected GA GNSS network were down¬ 
loaded from the GA archive 3 for 3 rd , 4 th , 6 th and 16 th 
March 2014 corresponding to Day of Year (DOY) 062, 
063, 065 and 075, respectively. The four days chosen 

3 ftp://ftp.ga.gov.au/geodesy-outgoing/gnss/data/ 
daily/yyyy /yyddd/xxxxdddO.yyd.Z The abbreviations yyyy 
and yy are the four and two digit year, ddd is the DOY, xxxx 
represents the four character GPS station id, d stands for 
Hatanaka compressed (Hatanaka, 2008) Receiver INdependent 
EXchange (RINEX) format (Gurtner &; Estery, 2007), and Z 
indicates compressed/zipped file. The Hatanaka compressed 
files can be decompressed by the software available at ftp:// 
terras.gsi.go.jp/software 


for this analysis were the first four nights of GLEAM 
observations for which suitable MWA data was avail¬ 
able. Also, by choosing data from the year 2014, data 
from the recently active GPS receiver MROl could be 
included. Only YAR3 is included in CODE analysis, il¬ 
lustrating why our analysis is required to establish dense 
GPS networks near the MWA. Table 4 presents the sum¬ 
mary of ionospheric weather on the selected four days, 
the parameters, F10.7 solar flux, and Planetary Kp in¬ 
dices are presented. The data presented in Table 4 in¬ 
dicate quiet ionospheric conditions, which are ideal for 
testing the methods described in this paper. 



• GA stations • MWA MWA IPP 


Figure 3.: Selected GPS station locations from Geo¬ 
science Australia’s network (red), MWA location (blue) 
and MWA IPP (green) for the four MWA observation 
nights (DOY 062, 063, 065 and 075 marked by 1 to 4). 


3 CODE IONEX maps 

The daily ionosphere maps from the CODE are based 
on a global network of ~200 GPS/GNSS stations. The 
line of sight GPS ionospheric delay is mapped to VTEC 
using a Modified Single Layer Model (MSLM) mapping 
function approximating the Jet Propulsion Laboratory 
(JPL) Extended Slab Model (ESM) (for ESM mapping 
function see, Coster et al., 1992). The MSLM 4 is given 
as follows: 

Fmslm = sin z' = ( R ■ sin(az)) . (14) 

cos z' \K+H ) 

where H = 506.7 km and a = 0.9782. 


4 www.aiub.unibe.ch/download/users/schaer/igsiono/doc/ 
mslm.pdf 
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Table 2: Kalman Filter parameters and definition 


Parameter 

Definition 

Units 

el 

Elevation angle of satellite 

(radians) 

t 

Time step 


tmax 

Maximum number of time steps 


Xo\0 

Initial estimate of unknowns 


-*^010 

Variance matrix of x 0 |o 


A T 

A 0|0 

Design matrix 


Qy 

Stochastic model of observations 


Uo\o 

Observation vector 



Transition matrix 



Predicted unknowns 


Xt-l\t-l 

Unknowns from previous time step 


s t 

Variance matrix of the system noise 


Vt 

Predicted residuals 


Qv t 

Variance matrix of predicted residual 


K t 

Kalman gain matrix 


x t \t 

Updated unknowns 


r x t \t 

Variance matrix of x t \t 



Table 3: Description of the selected GA network GPS/GNSS stations and the MWA 


Station Receiver type 


MROl TRIMBLE NETR9 
MTMA LEICA GRX1200+GNSS 
YAR3 LEICA GRX1200GGPRO 
WILU LEICA GRX1200+GNSS 
MWA 


Antenna type Observables Location Observing sessions 

used (degrees) (Year,Day of Year) 

TRM59800.00 LI, L2, Cl, P2 26.70° S 116.37° E 2014, 062, 063, 065, 075 

LEIAR25.R3 LI, L2, Cl, P2 28.11° S 117.84° E 2014, 062, 063, 065, 075 

LEIAR25 LI, L2, Cl, P2 29.04° S 115.34° E 2014, 062 a , 063, 065, 075 

LEIAR25.R3 LI, L2, Cl, P2 26.62° S 120.21° E 2014, 062, 063, 065, 075 

26.70° S 116.67° E 2014, 062, 063, 065, 075 


a Partial data available, from 00:00:00 UTC to 18:07:00 UTC 


Table 4: Daily solar and geomagnetic indices for the selected MWA observation days 


Year, DOY 

Solar Flux at 10.7 cm (F10.7“) 
10~ 22 Wm~ 2 H z~ L 

Planetary K index (Kp b ) 

3 hourly, from 00 to 24 UTC, ranging from 0-9 (low-high) 

2014, 062 

161 

3 

2 

2 

2 

2 

1 

1 1 

2014, 063 

158 

2 

1 

0 

2 

3 

2 

2 3 

2014, 065 

149 

1 

1 

1 

2 

2 

2 

3 1 

2014, 075 

136 

1 

1 

1 

0 

0 

0 

1 0 


a National Geophysical Data Center (NGDC), National Oceanic and Atmospheric Administration (NOAA) ftp://ftp.ngdc.noaa. 
gov/STP/GEOMAGNETIC_DATA/INDICES/KP_AP/2014 

b Space Weather Prediction Center (SWPC), NOAA, BOULDER, USA ftp://ftp.swpc.noaa.gov/pub/warehouse/2014 


The VTEC is modelled using spherical harmonic co¬ 
efficients of order and degree 15 in the solar magnetic 
reference frame as snapshots with an ionosphere re¬ 
freshing interval of 2 hours (Schaer et ah, 1996). The 
spatial resolution of CODE maps is 5°/2.5° in lon¬ 
gitude/latitude, respectively. The spherical harmonic 
model used by CODE to interpret the global ionosphere 
is described in Schaer et al. (1996). 


CODE maps are available as daily final solutions in 
CODE’S online archive 5 . Figure 4 presents the 2 hour 
snapshots for the day 03-03-2014 (DOY 062), which is 
one day of interest for the MWA observations. The av¬ 
erage CODE VTEC RMS is shown in Figure 5. The 

5 f tp: //f tp.unibe. ch/aiub/CODE /yyyy/CODGddd.yyl, where 
yyyy and yy are the four and two digit year, ddd is the DOY. 
The ionosphere maps are exchanged in IONosphere Map Ex¬ 
change Format (IONEX), see Schaer et al. (1998) for detailed 
description of IONEX data format and its interpolation. 
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Figure 4.: Global TEC (TECU) from CODE IONEX maps for DOY 062, year 2014. 


RMS of the VTEC fit is higher over oceans and regions 
with sparse GPS/GNSS receivers coverage, due to lim¬ 
ited data points available for the fit. 


4 MEASUREMENT OF THE 
IONOSPHERE USING MWA 
OBSERVATIONS 

4.1 Observations 

For comparison with our GPS modelling of the iono¬ 
sphere, we used observations from the GLEAM survey 
(Wayth et al. (2015); Hurley-Walker et al. in prep.). In 
this survey, the MWA observes in meridian drift-scan 
mode, where the telescope remains pointing at a single 
point on the meridian throughout the night. Four nights 
from March 2014 were chosen, when the telescope was 
pointed close to the zenith. The ionosphere pierce point 
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corresponding to the pointing centre of the telescope is 
shown for each night in (Figure 6). 

During GLEAM observations, the instrument cy¬ 
cles five frequency bands, centred on approximately 
(88MHz, 118MHz, 154MHz 185MHz and 215MHz), 
with a dwell time of 2 minutes on each band, each band 
having an instantaneous bandwith of 30.72 MHz. Each 
two-minute observation is then imaged with separate 
images being generated for 4 subbands (each having a 
bandwidth of 7.68 MHz). The position of a prepared list 
of bright sources was then determined for each of these 
images. 

4.2 Ionospheric Modelling 

By comparing the position of each source in all four sub¬ 
bands, it is possible to quantify the effect of the iono¬ 
sphere, separating it from instrumental and calibration 
effects. By making an least-squares fit to all four points, 
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Figure 5.: Average RMS (1-u uncertainties) in TECU, 
of CODE IONEX maps for DOY 062, year 2014, marked 
in white are the GPS/GNSS stations considered for the 
solution. 


A comprehensive analysis of this dataset using MWA 
observations is underway (Morgan et al. in prepara¬ 
tion). For this analysis, ionospheric gradient were esti¬ 
mated over all sources detected in each snapshot. This 
is shown in Figure 7 shows the gradient of this fit for 
each source, scaled by A 2 for the highest frequency off¬ 
set shown in the left panel. The fact that the average 
reduced y 2 for each observation is ~ 1, and the fact that 
A = 0 position of the sources remains at zero throughout 
the night, both serve to reinforce the hypothesis that the 
shift in sources is largely due to the ionosphere. 

For simplicity, only the lowest-frequency data were 
used, since these observations are the most sensitive to 
ionospheric effects, and still yield sufficient time resolu¬ 
tion. This yields two time series (on in the NS direction 
and one in the EW direction) representing the average 
shift due to the ionosphere of ^100 sources within a 
35°radius of the pointing centre. In order that these 
shifts could be compared directly with GPS measure¬ 
ments, both measurements were scaled to a common 
reference frequency of 150 MHz and the (angular) offset 
was converted to radians. 



Figure 6.: GA station locations (red) with satellite IPPs 
in earth fixed reference frame over a period of 24 hours 
(grey) and for 1 hour during MWA observations (black). 
MWA is marked in blue, whereas the IPP of MWA is 
shown in green. 


the contribution of the ionosphere (the gradient) can 
be determined for each source. The offset in position at 
each frequency is shown in Figure 8 for a single strong 
source. It can be seen that the change in apparent po¬ 
sition of the source depends precisely on A 2 , exactly as 
would be expected from ionospheric refraction. 
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5 RESULTS AND DISCUSSION 

5.1 Comparison of VTEC with CODE 
IONEX 

The VTEC was estimated by our software at a time 
resolution of 10 minutes as described in section 2 with 
software developed in Matlab for the location of the four 
GPS stations, MROl, MTMA, YAR3, and WILU. The 
CODE values of VTEC are available at intervals of 2 
hours. Values of the VTEC corresponding to our time 
resolution were interpolated from CODE VTEC maps 
(see Schaer et ah, 1998) for each of the four GPS lo¬ 
cations. Our estimated values of VTEC along with the 
CODE VTEC values are presented in Figure 9 for DOY 
062, year 2014. Figures 9(a), 9(b), 9(c), 9(d), present 
VTEC for MROl, MTMA, YAR3, and WILU respec¬ 
tively for DOY 062. Figure 10 shows the difference be¬ 
tween our VTEC estimates with respect to CODE. The 
F10.7 solar flux, presented in Table 4, is highest on DOY 
062 and lowest on 075. This is reflected in the VTEC 
values, they reach a maximum of ^68 TECU on DOY 
062 and ~57 TECU on DOY 075 for MROl, refer Fig¬ 
ures 9(a) and B.3(a). The 1 <t uncertainties in CODE 
maps reach a maximum of 8 TECU (Figure 5). The dif¬ 
ferences between CODE and our VTEC are found to he 
within the errors, with the differences ranging between 
-6 to 6 TECU for four different days of observations, 
(Figure 10). 
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Figure 7.: Offset (extrapolated to A = 0), gradient, and reduced \ 2 of a fit of source position offset as a function 
of A 2 . Each point is for a single observation, all quantities are averaged over all ( 100) sources detected in that 
observation. Left panels are for the East-West position offset (Right Ascension) right panels are for the North-South 
position offset (Declination). Red, Yellow, Green, Blue and Purple are for the 88 MHz, 118 MHz, 154MHz 185 MHz 
and 215 MHz bands respectively. Note that the gradients (in arcmin nr 1 ) have been multiplied by 4, representing 
an offset at a wavelength of 2 m (=150MHz). 


5.2 Comparison of Receiver DCBs 

The receiver DCBs given by c • d Tt 21 in equation (5) 
are the inter-frequency biases on code GPS data. Es¬ 
timation of receiver DCBs is important for correct esti¬ 
mation of ionospheric parameters from GPS/GNSS ob¬ 
servables (Gaposchkin & Coster, 1992; Sardon et al., 
1994; Teunissen & Kleusberg, 1998). The receiver DCBs 
were estimated as described in section 2. In addition, 
the Bernese GNSS data processing software (Beutler 
et al., 2007) was used to estimate receiver DCBs for 
the selected GA GPS/GNSS stations given in Table 3. 
Bernese 5.0 Precise Point Positioning (PPP) process¬ 
ing estimates one set of station specific ionosphere pa¬ 
rameters for the entire session as well as receiver DCBs 
using the script PPP_I0N (Beutler et al., 2007). With 
Bernese PPP it is possible to obtain centimetre level 
station positions using precise satellite orbits, satellite 
clock and Earth Orientation Parameters (EOP) along 
with the receiver DCB. 

The station specific ionosphere parameters are es¬ 
timated by Bernese PPP by forming frequency- 
differenced geometry-free observables. In this work, 
frequency-differenced observables were used to esti¬ 


mate station specific parameters. Unlike Bernese PPP 
where quotidian ionosphere parameters are estimated 
for a given session, we make use of a Kalman filter to 
estimate the ionosphere every 10 minutes and a single 
value of receiver-satellite DCB for an entire session. 

Since there exists a rank deficiency between the re¬ 
ceiver and satellite DCBs. To overcome this rank defi¬ 
ciency, CODE 6 assumes a zero-mean condition over all 

m 

the satellite DCBs, ]T) DCB S = 0. By assuming a zero- 

mean condition it implies that the DCB results may be 
shifted by a common offset value (see Beutler et al., 
2007, Chapter 13), which is a function of total number 
of satellites m considered for the solution. Hence an in¬ 
dependent 5-basis was formed for estimation of DCBs. 

Our estimated receiver DCBs and Bernese PPP es¬ 
timates are compared in this section. CODE also pro¬ 
vides estimates of receiver DCBs which are estimated 
while performing a global fit for ionosphere parameters. 
However, of the selected stations used in this study, 
only receiver DCBs for station YAR3 are available from 
CODE. 


6 http://igscb.jpl.nasa.gov/igscb/center/analysis/ 
archive/code_20080528. acn 
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Figure 8.: Source position offset against A 2 (m 2 ) for a strong source. Each line represents a measurement of the 
source in each of four subbands in a single 2-minute observation. There are many lines since the source in multiple 
observations as it passes through the field of view. Left panels are for the East-West position offset (Right Ascension) 
right panels are for the North-South position offset (Declination). Red, Yellow, Green, Blue and Purple are for 
the 88 MHz, 118 MHz, 154 MHz 185 MHz and 215 MHz bands respectively. More significant detections are given a 
darker colour. 


Figure 11 shows our estimated receiver DCBs (blue), 
with Bernese (red) and from CODE (cyan) for DOY 
062. Table 5 presents the differences between our es¬ 
timated receiver DCBs and Bernese and CODE val¬ 
ues. In comparison with Bernese estimated DCBs, dif¬ 
ferences are between -0.475 to -1.311 ns which corre¬ 
sponds to -1.356 to -3.743 TECU. Whereas comparing 
to CODE DCBs, the DCB for YAR3 differed by 0.055 
ns to 0.553 ns, 0.157 to 1.579 TECU. 

Hong et al. (2008) estimated the receiver DCBs by 
initially estimating single differenced DCBs. Further, 
by finding the time to for which the single difference 
geometric range is zero, absolute receiver DCBs were 
computed. Hong et al. compared the estimated receiver 
DCBs with Bernese values and found the maximum 


difference to lie around 15cm or 0.5ns. Arikan et al. 
(2008) estimated the receiver DCBs using IONOLAB- 
BIAS method and compared them with CODE esti¬ 
mates. Arikan et al. found the differences in receiver 
DCBs to lie between -0.552ns and 0.110ns for different 
receivers. 


5.3 Comparison of GPS ionosphere gradients 
with MWA observations 

The gradients in the East-West (EW) and North-South 
(NS) directions were computed from the ionosphere first 
order coefficients given in equation (11) for each of the 
selected stations and compared with MWA observed 
gradients. Figure 12 presents each of the EW and NS 
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MROl MTMA 



(a) VTEC at station MROl (b) VTEC at station MTMA 


YAR3 WILU 



(c) VTEC at station YAR3 (d) VTEC at station WILU 


Figure 9.: VTEC at stations MROl, MTMA, YAR3, and WILU estimated using the method described in the text 
(blue curve) and CODE IONEX (red curve) on DOY 062, year 2014. 


Table 5: Differences between our estimated receiver DCBs with Bernese and CODE for DOY 062, 063, 065 and 
075, year 2014. 


Station Difference (ns) Difference (ns) 

w.r.t Bernese w.r.t CODE 

Year 2014, DOY Year 2014, DOY 



062 

063 

065 

075 

062 

063 

065 

075 

MROl 

-1.311 

-0.861 

-0.651 

-0.648 

- 

- 

- 

- 

MTMA 

-0.640 

-0.760 

-0.629 

-0.882 

- 

- 

- 

- 

YAR3 

-1.262 

-0.475 

-0.778 

-0.756 

0.553 

0.230 

0.095 

0.055 

WILU 

-0.924 

-0.806 

-0.711 

-0.980 

- 

- 

- 

- 


gradients for all the GA GPS stations and MWA. Fig¬ 
ures 12(a), 12(c), 12(e), and 12(g), show the EW gra¬ 
dients for DOY 062, 063, 065 and 075 of year 2014, 
respectively. The NS gradients are presented in Figures 
12(b), 12(d), 12(f), and 12(h) for DOY 062, 063, 065 and 
075 of year 2014, respectively. In each of the subplots, 
along with the GPS ionospheric gradients, the MWA 
observed gradients are shown. 


Table 6 presents the correlation between GPS and 
MWA gradients in the EW and NS directions for the 
four days of observations. The IPP separations in lon¬ 
gitude (|AAjpp|) and latitude (|A</>/pp|) for each of 
the GPS stations and MWA are presented in Table 6. 
The correlation between the GPS and MWA ionosphere 
gradients is computed using Pearson’s coefficient of cor- 
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(a) Differences in VTEC on DOY 062 (b) Differences in VTEC on DOY 063 



(c) Differences in VTEC on DOY 065 
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Figure 10.: Differences in VTEC w.r.t CODE for stations MROl, MTMA, YAR3, and WILU for DOY 062, 063, 
065 and 075 year 2014. 


relation, r, for an assumed mean method along with the 
standard error, ay, refer Fisher (1936). 

There is a high correlation in the EW and NS gradi¬ 
ent between GPS and MWA for most of the GPS sta¬ 
tions for all the days (Table 6). The NS gradients had a 
weak correlation with YAR3 for most days. The corre¬ 
lation was highest on DOY 075 while the IPP was clos¬ 
est among the four days of observation (Table 6). The 
EW gradient showed consistent good correlation with 
MROl GPS stations, whereas correlation with WILU 
seemed to be most inconsistent (Table 6). A general 
trend seemed to show that the EW gradient was propor¬ 
tional to the longitudinal difference and the NS gradient 
to the latitude difference. For DOY 075 while the solar 
activity was the lowest among the selected days (Table 
4), the EW gradient was found to be high for all the 
stations. NS gradient did not have a similar behaviour. 

Table 7 and Figure 12 present the comparison of 
zenith EW and NS gradients between GPS stations. 
Correlation for the EW and NS gradients, presented 
in Table 7, were computed between each of the selected 


GPS stations for the time window of MWA observations 
(marked by the red line in Figure 12). Table 7 sum¬ 
marises the inter-station EW and NS gradient correla¬ 
tions. Inter-station correlation is strong for the EW and 
NS gradient for almost all the days between all the sta¬ 
tions. The EW gradient is consistently strong between 
all the stations, however the NS gradient is found to be 
weakest between WILU and YAR3 which have a distin¬ 
guished latitudinal and longitudinal separation (Figure 
3). 


6 CONCLUSIONS 

This work presents a single-station approach to esti¬ 
mate VTEC ionosphere gradients. The VTEC and 
ionosphere gradients were estimated at intervals of 10 
minutes. The ionosphere gradients in the EW and NS 
direction at the GPS station locations are in good agree¬ 
ment with the MWA observed gradients. 

The ionosphere gradient analysis presented in this re¬ 
search brings forth various questions, namely, the vari- 
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Table 6: Correlation between the GPS and MWA observed gradients in EW {tew) and NS ( tns ) components, 
its standard error (ay), and IPP separations in longitude (|AA/pp|) and latitude (\Acf>ipp\) of GPS stations and 
MWA (AIPP) for DOY 062, 063, 065 and 075, year 2014. 

Station r EW (oy EW ) |AA/pp| (degrees) 

Year 2014, DOY Year 2014, DOY 



062 

063 

065 

075 

062 

063 

065 

075 

MROl 

0.79(0.05) 

0.73(0.06) 

0.66(0.08) 

0.93(0.02) 

0.03 

0.03 

0.03 

0.03 

MTMA 

0.72(0.07) 

0.65(0.08) 

0.71(0.07) 

0.94(0.02) 

1.17 

1.17 

1.17 

1.17 

YAR3 

0.83 a (0.05) 

0.73(0.06) 

0.77(0.06) 

0.92(0.02) 

1.32 

1.32 

1.32 

1.32 

WILU 

0.60(0.09) 

0.83(0.04) 

0.54(0.10) 

0.94(0.01) 

3.54 

3.54 

3.54 

3.54 



tns (oVjvs) 


\Ac/) I pp\ (degrees) 

MROl 

0.84(0.04) 

0.78(0.05) 

0.87(0.03) 

0.69(0.07) 

0.10 

0.65 

1.51 

0.65 

MTMA 

0.77(0.06) 

0.75(0.06) 

0.74(0.06) 

0.46(0.11) 

1.31 

2.07 

2.93 

0.77 

YAR3 

0.21 a (0.16) 

0.45(0.11) 

0.50(0.10) 

0.35(0.12) 

2.24 

3.00 

3.86 

1.70 

WILU 

0.89(0.03) 

0.77(0.05) 

0.86(0.04) 

0.64(0.08) 

0.17 

0.58 

1.43 

0.72 


“Partial data available, from 00:00:00 to 18:07:00 UTC (MWA observation window = ~ 11 to 21 UTC) 

Table 7: Inter-station correlation for the EW and NS gradients (r), its standard error (ay), and IPP separations in 
longitude (|AA/pp|) and latitude (|A0/pp|) between GPS stations 


Station 


tew 

(°Vew ) 


AA/pp| (degrees) 

062 

Year 2014, DOY 

063 065 

075 


MROl-MTMA 

0.88(0.03) 

0.97(0.01) 

0.94(0.02) 

0.99(0.002) 

1.21 

YAR3-MTMA 

- 

0.83(0.04) 

0.94(0.01) 

0.97(0.01) 

2.50 

MTMA-WILU 

0.90(0.03) 

0.91(0.02) 

0.92(0.02) 

0.98(0.01) 

2.37 

MR01-YAR3 

- 

0.89(0.03) 

0.91(0.02) 

0.97(0.01) 

1.29 

MROl-WILU 

0.88(0.03) 

0.96(0.01) 

0.93(0.02) 

0.97(0.01) 

3.58 

WILU-YAR3 

- 

0.86(0.04) 

0.84(0.04) 

0.94(0.02) 

4.87 



tns 

(°V; \is) 


\A4>ipp\ (degrees) 

MROl-MTMA 

0.94(0.02) 

0.92(0.02) 

0.95(0.01) 

0.89(0.03) 

1.42 

YAR3-MTMA 

- 

0.79(0.05) 

0.88(0.03) 

0.93(0.02) 

0.93 

MTMA-WILU 

0.81(0.05) 

0.87(0.03) 

0.95(0.01) 

0.90(0.03) 

1.49 

MR01-YAR3 

- 

0.71(0.07) 

0.83(0.04) 

0.84(0.04) 

2.35 

MROl-WILU 

0.84(0.04) 

0.91(0.02) 

0.96(0.01) 

0.96(0.01) 

0.07 

WILU-YAR3 

- 

0.53(0.10) 

0.74(0.06) 

0.79(0.05) 

2.42 


PASA (2015) 

doi: 10.1017/pas.2015.xxx 






















MWA ionosphere calibration 


17 


2 5 r 


20 - 


15 - 


10 - 


25 


20 - 


15 - 


10 - 


■■ 



^:jps Manmnl Icons \ 

MROl 

WILU YAR3 MTMA 

Receiver name 

(a) Receiver DCBs for DOY 062, year 2014 


I' 


WILU YAR3 

Receiver name 


25 


20 


15 


10 


25 


20 


15 


10 


1 ' 




^GPS MBernesi ZlCODE ( 

MROl WILU YAR3 MTMA 

Receiver name 

(b) Receiver DCBs for DOY 063, year 

2014 



' 


I GPS ^MBernesF^lCODE [. 


!■ 


I GPS MWBernesj ICODE [ 


WILU YAR3 

Receiver name 


(c) Receiver DCBs for DOY 065, year 2014 


(d) Receiver DCBs for DOY 075, year 2014 


Figure 11.: Comparison of our estimated receiver DCBs (blue) with Bernese estimates (red) and CODE DCBs 
(cyan) for DOY 062, 063, 065 and 075, year 2014. 


ation of EW and NS gradients with respect to IPP 
separation, and the elevation dependency reflected by 
the correlation of gradients between the GPS station, 
MROl, and MWA. With the limited data set, these 
questions cannot not be resolved comprehensively. The 
future work will focus on including more MWA observa¬ 
tion and will attempt to closely probe the above ques¬ 
tions in a statistical sense. 

With our single-station approach, the GPS receiver 
DCBs can be accurately estimated for any available 
GPS receiver. Thus our method can be applied to local 
GPS receiver data for which the DCBs are not publicly 
available. 

To develop a regional model for the ionosphere, a 
multi-station approach needs to be adopted. Further¬ 
more, the spatial resolution of ionosphere gradients is 
closely related to the scale of the GPS/GNSS receiver 


network on the ground. The scales at which MWA sees 
the ionosphere lie between 10 - 100 km. In order to es¬ 
timate the ionosphere gradients on such scales, dense 
GPS networks of the order of the ionosphere scales seen 
by the MWA need to be present around MWA site. The 
existing GPS network near the MWA site is of scale 
150 - 250 km, which limits the GPS based ionosphere 
research possible for the MWA. 

In addition to GPS satellite data, other satellite sys¬ 
tems can be used to density the IPPs over the study 
area. Other global navigation systems like GLONASS 
(GLObal NAvigation Satellite System), launched by 
Russia, BeiDou from China, and the European Union’s 
Galileo could be adopted in our model. Of these three 
global navigation systems, only the GLONASS system 
is close to completely operational, it currently has 23 
operational satellites in orbit. GA receivers are able 
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Hour (UTC) 

(a) EW gradient for DOY 062 


Hour (UTC) 

(b) NS gradient for DOY 062 




(c) EW gradient for DOY 063 


(d) NS gradient for DOY 063 



Hour (UTC) 

(e) EW gradient for DOY 065 


Hour (UTC) 

(f) NS gradient for DOY 065 




(g) EW gradient for DOY 075 


(h) NS gradient for DOY 075 
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to capture data from the GLONASS system, and us¬ 
ing GLONASS data along with the GPS data has been 
shown to improve ionospheric modelling (Coster et al., 
1999). 

Future work will focus on developing a regional iono¬ 
sphere model using data from both GPS and GLONASS 
satellite systems. 

We find that the Australian SKA site (where the 
MWA is located) is well suited for low-frequency as¬ 
tronomy. Appendix A shows that the conclusions drawn 
in Sotomayor-Beltran et al. (2013) regarding the suit¬ 
ability of MRO for low-frequency radio astronomy are 
a result of an incorrect interpretation of CODE maps. 
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A CODE IONOSPHERE MAP 

INTERPRETATION DISCREPANCY 

The CODE maps presented in Sotomayor-Beltran et al. 
show a behaviour of the equatorial anomaly different to 
what is seen in this work. To underline this inconsistency, 
the CODE maps for the day 11 April 2011 presented by 
Sotomayor-Beltran et al. (Figure Ala) were plotted with 
our software, Figure Alb. The possible source of inconsis¬ 
tency lies in the way CODE maps are plotted by Sotomayor- 
Beltran et al.. On close inspection it is clear that the lat¬ 


itudes provided in CODE IONEX files are inverted in the 
plot presented by Sotomayor-Beltran et al.. This results in 
the equatorial anomaly appearing to pass directly over the 
MRO, which is not the case. This conclusion is supported in 
the literature, Seeber (2003, Figure 7.52), Kennewell et al. 
(2005). 


06 UT |hr| 



(a) CODE IONEX plot in Sotomayor-Beltran 
et al. (2013) 


0 6 UTC 



(b) CODE IONEX plot by our software 


Figure Al.: CODE IONEX plot for 06 UTC by (a) 
Sotomayor-Beltran et al. (2013) and (b) by our soft¬ 
ware for 11 April 2011 
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Figure B.I.: VTEC at stations MROl, MTMA, YAR3, and WILU estimated using the method described in the 
text (blue curve) and CODE IONEX (red curve) on DOY 063, year 2014. 
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Figure B.2.: VTEC at stations MROl, MTMA, YAR3, and WILU estimated using the method described in the 
text (blue curve) and CODE IONEX (red curve) on DOY 065, year 2014. 
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Figure B.3.: VTEC at stations MROl, MTMA, YAR3, and WILU estimated using the method described in the 
text (blue curve) and CODE IONEX (red curve) on DOY 075, year 2014. 
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